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Abstract 

A key problem of eukaryotic cell motility is the signaling mechanism of chemoattractant gradient sens- 
ing. Recent experiments have revealed the molecular correlate of gradient sensing: Frontness molecules, 
such as PI3P and Rac, localize at the front end of the cell, and backness molecules, such as Rho and 
myosin II, accumulate at the back of the cell. Importantly, this frontness-backness polarization occurs 
"spontaneously" even if the cells are exposed to uniform chemoattractant profiles. The spontaneous po- 
larization suggests that the gradient sensing machinery undergoes a Turing bifurcation. This has led to 
several classical activator-inhibitor and activator-substrate models which identify the frontness molecules 
with the activator. Conspicuously absent from these models is any accounting of the backness molecules. 
This stands in sharp contrast to experiments which show that the backness pathways inhibit the frontness 
pathways. Here, we formulate a model based on the mutually inhibitory interaction between the front- 
ness and backness pathways. The model builds upon the mutual inhibition model proposed by Bourne 
and coworkers (Xu et al, Cell, 114, 201-214, 2003). We show that mutual inhibition alone, without 
the help of any positive feedback, can trigger spontaneous polarization of the frontness and backness 
pathways. The spatial distribution of the frontness and backness molecules in response to inhbition and 
activation of the frontness and backness pathways are consistent with those observed in experiments. 
Furthermore, depending on the parameter values, the model yields spatial distributions corresponding to 
chemoattraction (frontness pathways in-phase with the external gradient) and chemorepulsion (frontness 
pathways out-of-phase with the external gradient). Analysis of the model suggests a mechanism for the 
chemorepulsion-to-chemoattraction transition observed in neurons. 

Keywords: Eukaryotic cells, chemotaxis, gradient sensing, directional sensing, spontaneous polarization, 
Turing instability 

1 Introduction 

When motile cells are exposed to a chemoattractant gradient, they develop a morphological polarity con- 
sisting of a distinct front and back p]. The formation of the morphological polarity is driven by the spatial 
segregation of distinct sets of intracellular molecules to the front and the rear of the cell (see Figure QJi and 
refs. ElElllEl). The frontness molecules, which include Cdc42, Rac, PI3P, PI3K, Arp2/3, and F-actin, 
localize to the front of the cell where they coordinate the extension of an actin-rich extension. The backness 
molecules, which include Rho, Rho kinase, PTEN, and myosin II, migrate to the rear of the cell where they 
are thought to activate cell contraction. 

In general, the spatial segregation of the frontness/backness molecules and the resultant morphological 
polarization occurs even if the cells are exposed to a uniform chemoattractant profile (see Figure QJ) and 
OdIHI). This phenomenon has been called spontaneous polarization to emphasize the fact that the cells 
polarize despite the absence of a perceptible external cue 

The existence of spontaneous polarization is reminiscent of the Turing instability in reaction-diffusion 
systems pHl- Thus, it has led to several models that view spontaneous polarization, either explicitly or 
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Figure 1: Spatial segregation of frontness and backness molecules, (a) When neutrophils are exposed to 
a gradient of fMLP, Rho (represented by its marker, RBD) and PTEN colocaHze at the back of the cells, 
and Cdc42 (represented by its marker, CBD) localizes to the front of the cells (from |5J). (b) The spatial 
segregation occurs even in the absence of a chemoattractant gradient. When neutrophil- like HL-60 cells are 
exposed to a uniform concentration (100 nM) of fMLP, actin polymers (shown in red) localize at the front 
of the cell, and Rho (shown in green) locaHzes at the back of the cell (from ^). 



implicitly, as the onset of a Turing bifurcation ^l] El El El El According to these models, the cell is 
in a stable homogeneous steady state in the absence of chemoattractant. However, exposure of the cell to a 
sufficiently large uniform chemoattractant concentration pushes it past a Turing bifurcation point, where the 
homogeneous steady state is unstable with respect to certain non-homogeneous perturbations. The inevitable 
presence of noise is then sufficient to drive the cell towards a nonhomogeneous steady state corresponding 
to the polarized state of the cell. 

As attractive as these models of spontaneous polarization may be, there is a significant gap between 
the theory and experiments. Indeed, all the models consist of one or more slow-diffusing activators whose 
synthesis is autocatalytic, and which in consequence, tend to grow and spread across the entire cell. The 
unconstrained growth and dispersion of the activator is restricted by hypothesizing the existence of a diffusible 
inhibitor (which is a by-product of activator synthesis and impedes the growth of the activator) or substrate 
(which is consumed during activator synthesis and stimulates the growth of the activator). The diffusible 
activator/substrate ensures that the growth of the activator(s) remains confined to a localized region of the 
cell membrane. This region is identified with the front of the cell, and the rest of the cell membrane, suffering 
from an activator deficit, is presumed to constitute the back of the cell. The predictions of these models 
are partially consistent with experiments involving the activation or inhibition of the frontness molecules. 
Indeed, if PI3K is inhibited, the gradient of the frontness molecules such as PI3P progressively decreases, 
until at sufficiently high levels of inhibition, there is no gradient at all (Figure |2t)- On the other hand, 
when Rac is overexpressed, high levels of frontness molecules, PI3P and Rac, are found all over the cell 
membrane (Figure (Jb). Both these features are reproduced by the activator- inhibitor class of models |12II14I 
El- However, these models cannot explain the spatial distribution of the backness molecules in response to 
inhibition of the frontness molecules. Specifically, if the activity of Cdc42 is suppressed, Rho is found not only 
at the back but also at the front of the cell (FigureEt)- Likewise, if the Gi proteins are inhibited with pertussis 
toxin (PTX), a uropod-like structure forms at the up-gradient edge of the cell (Figure Oa). The models also 
offer no insight into experiments involving activation or inhibition of the backness molecules. For instance, 

^Alternative models of gradient sensing in which spontaneous polarization plays no role have also been proposed 1161 1171 
1181 1191 12?1] . For the most part, they are motivated by experimental systems such as PDGF-stimulated fibroblasts |2n| and 
latrunculin-treated Dictyostelium cells 1161 . which do not display spontaneous polarization. In these models, polarization 
occurs only in the presence of external gradients. 
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Figure 2: Spatial redistribution of the frontness pathways in response to inhibition and activation of the 
frontness pathways, (a) When HL-60 cells are pretreated with intermediate concentrations (100-200 /iM) 
of LY294002, an inhibitor of PI3K, before their exposure to fMLP, the gradient of PH-Akt, a marker for 
PI3P, decreases. The figures on the left and right show the spatial distribution and fluorescence intensity of 
PH-Akt in control and LY294002-treated cells, respectively. At high concentrations (300 /xM) of LY294002, 
PH-Akt does not polarize at all (from 0). (b) When HL-60 cells transfected with constitutively active 
Rac are exposed to fMLP, the frontness molecules, PI3P (top) and active Rac (bottom), spread all over the 
membrane (from ['21\). 



when Rho is activated, PI3P fails to polarize (Figure 0Ji). Conversely, when Rho is inhibited, PI3P spreads 
all over the membrane, and the cell extends a single broad pseudopod or multiple pseudopods (Figure ^b)- 
These results are beyond the scope of the activator-inhibitor and activator-substrate models because the 
backness molecules are not even acknowledged as legitimate variables — they are implicitly assumed to 
somehow settle down in regions uninhabited by the frontness molecules. This stands in sharp contrast to the 
experimental data which shows that backness pathways downregulate the frontness pathways: Inhibition of 
Rho kinase increases Rac activity 2-3 fold 

One hypothesis regarding the interaction between the frontness and backness pathways is that they inhibit 
each other jj. Based on extensive experiments with neutrophils and neutrophil-like HL-60 cells |22l l6l 171 l2Tl 
15, Bourne and coworkers arrived at the kinetic scheme shown in Figure Et, which they describe as follows 

Briefly, the attractant binds to a G protein-coupled receptor (R) , which in turn activates different 
trimeric G proteins to generate two divergent, opposing signaling pathways, which promote polar- 
ized frontness and backness, respectively. In the frontness pathway, Gi, PI3Ps, and Rac promote 
de novo formation of actin polymers. One or more positive feedback loops in this first pathway 
mediate locaHzed increases in sensitivity to attractant: one of these requires polymerized actin, 
while Rac or Cdc 42 may in addition enhance PI3P accumulation more directly, via an alterna- 
tive pathway (dotted curved line in Figure 7 [reproduced here as Figure]^]). Backness signals, 
generated by G12 and G13, depend on activation of a Rho-dependent pathway that stimulates 
activation of myosin II, formation of contractile actin-myosin complexes, and myosin-dependent 
inhibition of Rac- and PI3P-dependent responses. Backness signals inhibit frontness signals, and 
vice versa (dashed straight Hnes in Figure 7). 

They go on to explain spontaneous polarization in terms of this kinetic scheme as follows 

The more or less symmetrically distributed actin ruffles and PI3P accumulation seen at early 
times (e.g., 30 s) after appHcation of a uniform stimulus presumably mask a fine-textured mosaic 
of interspersed backness and frontness signals, some triggering activation of PI3Ks, Rac, and 
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Figure 3: Spatial redistribution of the backness pathways in response to inhibition of the frontness pathways, 
(a) When neutrophils devoid of a-Pix, an activator of Cdc42, polarize spontaneously, the backness component 
RhoA is found at the back as well as the front of the cell (along with actin polymers) . The panels on the left 
and right show the distribution of F-actin (green) , RhoA (red) and the overlay in wild-type and a-Pix-null 
cells (from [5]). (b) HL-60 cells treated with PTX, a potent inhibitor of Gi proteins, extend a uropod-like 
structure at the edge exposed to the chemoattractant source (from fj). In control cells (top panel), a broad 
"head" is extended toward the fMLP source. In PTX-treated cells (middle and bottom panel), a knob-like 
structure, similar to the "tail" of control cells, is extended toward the fMLP source. 
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Figure 4: Spatial redistribution of the frontness and backness pathways in response to inhibition and ac- 
tivation of the backness pathways (from Q). (a) When HL-60 cells transfected with constitutively active 
Rho are exposed to fMLP, there is no polarization of PI3P. The top and lower panel show the evolution of 
the spatial distribution of PI3P in wild-type and transfected cells, (b) When HL-60 cells transfected with 
dominant negative Rho are exposed to a uniform concentration of fMLP, PH-Akt, a marker for PI3P, spreads 
all over the membrane, and multiple pseudopods develop. 



actin polymerization, others promoting activation of Rho and myosin. LocaHzed mechanical 
incompatibility of the two cytoskeletal responses, combined with the ability of each to damp 
signals that promote the other (dashed inhibitor lines in Figure 7) , then gradually drive them to 
separate into distinct domains of the membrane. 

The goal of this work is to formulate a mathematical model of spontaneous polarization based on the foregoing 
mechanism, namely, mutual inhibition of the frontness and backness pathways. 

It turns out that the Bourne model imposes two requirements that are incompatible with Turing insta- 
bilities in a two-component system consisting of frontness and backness pathways. Specifically, the mutual 
inhibition between the frontness and backness pathways must be sufficiently strong to ensure that the homo- 
geneous steady state is unstable in the presence of diffusion, and yet weak enough to guarantee its stability 
in the absence of diffusion. These two requirements cannot be satisfied simultaneously in a two-component 
system, a point that will be discussed in Sectional 

Thus, we are led to consider the modification of the Bourne model shown in Figure Eb. It contains two 
variables, J72, C^3, representing the frontness and backness pathways, respectively, whose activation is driven 
by receptor ligation, and which inhibit each other in a concentration-dependent manner. It is assumed fur- 
thermore that the mutual inhibition between U2 and U3 is so strong that they are mutually incompatible 
(in a sense that will be made mathematically precise below). We refer to these two variables as activators. 
The model differs from the Bourne picture inasmuch as it assumes the existence of a diffusible inhibitor, 
denoted Ui, whose synthesis is promoted by both activators, but which, in turn, inhibits both activators. 
The inhibitor serves to stabilize the coexistence of the otherwise incompatible activators: Transient increases 
in activator concentrations at localized "hotspots" are efficiently suppressed by the concominant increase 
in the concentration of the inhibitor. However, the mobility of the inhibitor imposes constraints upon its 
stabilizing effect because it tends to diffuse away from the "hotspot." This is not an issue at low chemoat- 
tractant concentrations, for under these conditions, the activator concentrations, and hence, their mutual 
incompatibility, is so small that despite the high diffusibility of the inhibitor, it successfully damps fiuctua- 
tions of the activator levels. However, at high chemoattractant concentrations, the activator concentrations 
and their mutual incompatibility are so large that the diffusible inhibitor fails to suppress the two activators 
sufficiently. The mutual incompatibility of the two activators now overcomes the mollifying effect of the 
inhibitor, and the activators segregate spatially into separate domains. 
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Figure 5: Model schemes, (a) Model proposed by Bourne and co-workers for spontaneous polarization of 
neutrophils and HL-60 cells (from 0]). The chemoattractant binds to G-protein-coupled receptors, R. The 
Gi-coupled receptors activate the frontness pathway which includes PI3P, Rac, and actin polymers. The 
G12/13-coupled receptors activate the backness pathway which includes Rho and myosin. The frontness 
pathway is subject to positive feedback from actin polymers (full line from actin polymers to PI3P) and 
possibly Rac (dashed line from Rac to PI3P). The frontness and backness pathways inhibit each other 
(dashed lines with a bar at the end), (b) The model scheme considered in this paper. Receptor activation 
stimulates the synthesis of frontness and backness subsystems, denoted U2 and U3, respectively. The frontness 
and backness subsystems stimulate their own synthesis (dashed lines), but inhibit each other both directly 
(line from U2 to C/3, and vice versa), and indirectly through production of the inhibitor, Ui. 
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The mathematical model in this paper quantifies the above physical argument. We show that the model 
yields spontaneous polarization at sufficiently high active receptor levels. Moreover, if the frontness or 
backness pathways are inhibited or activated, they redistribute spatially in a manner consistent with the 
experiments described above. Finally, we show the model displays steady states corresponding to both 
chemoattraction and chemorepulsion. Interest in this phenomenon is motivated by the fact that when 
neurons are exposed to activators of the cGMP pathway, chemorepulsion turns into chemorepulsion 23J. 
Based on the analysis of the model, we suggest that such transitions can be triggered by altering the balance 
of power in the mutually inhibitory interactions between the frontness and backness pathways. 

The mechanism of spontaneous polarization in this model is distinct from that in activator-inhibitor 
or activator-substrate models. UnHke these models, the spatial segregation of the two activators is driven 
entirely by their mutual inhibition — positive feedback plays no role. In reality, both positive feedback and 
mutual inhibition cooperate to produce symmetry-breaking (Figure Et)- The contribution of this work is to 
highlight the distinct role of mutual inhibition, a feature that emphasized in the experimental literature |24l 

but absent from previous activator- inhibitor and activator-substrate models. 

The paper is organized as follows. In Section [2l we define the model and simulate the experiments 
described above. In Sectional we elaborate on the physics underlying the model. Finally, we summarize the 
conclusions. 



2 Results 

2.1 The model 

We assume that 

1. Both activators promote the synthesis of the inhibitor, which in turn, degrades by a first-order process. 
Thus, the net rate of synthesis of J7i is 

-Riui + A12U2 + A13U3. 

2. The synthesis of U2 is receptor-mediated. It is autocatalytic at low concentrations and self-limiting at 
high concentrations, i.e., 

where r denotes the receptor activity, and the second term represents a self-limiting process that 
prevents synthesis rate from increasing beyond bounds. 

3. The synthesis of U2 is inhibited in a concentration-dependent manner by the common inhibitor, Ui, 
and the other activator, U3. Assuming that these interactions follow bimolecular kinetics, the net rate 
of synthesis of U2 is 

R2rU2 - A22u\ - A21U1U2 - A23U2U3. 

4. The synthesis of U3 follows kinetics similar to those of U2, i.e., its net rate of synthesis is 

i?3rM3 - A33U3 - A31U1U3 - A32U2U3 

Given these assumptions, we arrive at the equations 

dUl „ dUl r, . . 

= -Di^-^ - + A12U2 + A13U3 (1) 

~dT ^ ^ ^^'^^ ^ ^21"! - ^23^3 - A22U2) U2 (2) 

'dT ^ ^^'dX^ ~^ ^^^^ ~ ^3lUl - ^32^2 - A33U3) U3 (3) 

where X denotes the spatial coordinate, T denotes time, and ui, M2, U3 denote the concentrations of J7i, C/2, C/3, 
respectively. We assume the Neumann boundary conditions 

dui 

7— = at a; = 0, L 
oX 
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Figure 6: The bifurcation diagram for the homogeneous steady states (a) In the presence of the inhibitor 
and (b) in the absence of the inhibitor. The stable and unstable steady states are represented by full and 
open circles, respectively. 

since they imply no interaction with the environment, which is consistent with our desire to study spontaneous 
(autonomous) polarization. 

It is convenient to rescale the equations. If we define 



t 



aiiT,x = — , 



di 



di-i — 



Ri 



,P2 



P3 



Ri 

Ri 



we obtain the scaled equations 

du\ 
'dt 

dU2 

'dt 

dU3 

dt 

and boundary conditions 



, dui 



du2 



p2{r - a2lUl - 022^2 - a23'W3) W2 

Psir - asiui - 032^2 - 033^3) U3 



(4) 
(5) 
(6) 



dx 



at X = 0, 1. 



Our goal is to study the variation of the steady states as a function of the receptor activity, r. We shall 
show, in particular, that the model has a homogeneous steady state at which U2 and coexist, but it 
becomes Turing unstable at a sufficiently large r. The non-homogeneous steady state emerging from the 
Turing bifurcation is such that U2 and M3 are spatially segregated. 



2.2 The homogeneous steady states 

A Turing instability occurs when a non-homogeneous steady state bifurcates from a homogeneous steady 
state that is stable in the absence of diffusion. It is therefore useful to characterize the conditions for stability 
of the homogeneous steady states in the absence of diffusion. 
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The model has three types of homogeneous steady states: the trivial steady state, ui — U2 — = 0, 
denoted Eqqq', the semitrivial steady states, ui,M2 > 0,U3 = 0, and mi,U3 >,U2 = 0, denoted Euq and 
Eioi, respectively; and the nontrivial steady state, ui, 1*2,^3 > 0, denoted Em. The key results, depicted 
graphically in Figure Et, are as follows (see Appendix El for details): 

1. The trivial steady state, £'000, exists for all r > 0, and is always unstable. 

2. The semi-trivial steady state, Euq, which lies on the uiM2-plane, exists for all r > 0. It is stable if and 
only if 

/332^^>1. 

a22 

Here, 0:22 = 022 + 021012 is a measure of the extent to which U2 inhibits itself both directly (022) and 
indirectly through production of Ui (021012). Likewise, 032 = 032 + 031012 is a measure of the extent 
to which U2 inhibits U3 directly (032) and indirectly through production of Ui (031012). The foregoing 
stability condition then says that Euq is stable if and only if U2 inhibits U3 more than it inhibits itself. 

3. The semi-trivial steady state, i?ioi, which lies on the uiW3-plane, exists for all r > 0. It is stable if and 
only if 

P23 = > -L, 

^33 

where 0:23 = 023 -I- 021013,033 = 033 + 031013 are measures of the extent to which U3 inhibits U2 and 
itself, respectively, and (323 is measure of the cross-inhibition relative to the self-inhibition. 

4. The coexistence steady state, Em, exists if and only if 

f323,f332 < 1 or /323,/332 > 1- 

It is stable only if the mutual inhibition of U2 and C/3 is sufficiently weak, i.e., I323P32 < 1- 

Figure is reminiscent of the bifurcation diagram for the Lotka-Volterra model j^. This is not surprising 
because U2 and U3 obey Lotka-Volterra dynamics in the absence of the inhibitor. Indeed, the equations 
obtained by letting ui = in equations (j^El are identical to the Lotka-Volterra model for two competing 
species. It follows that the dynamics of U2 and U3 in the absence of the inhibitor are described by a 
bifurcation diagram very similar to Figure Et, the only difference being that (323 and (332 must now be 
replaced by 623 = 023/022, ^32 = 032/033 (see Figure Eb)- We shall appeal to this fact below. 



2.3 Turing instability of homogeneous steady states 

Since i^ooo is always unstable, it can never undergo a Turing bifurcation. However, the semi-trivial and 
non-trivial steady states are stable for all sufficiently small r > 0. The question then arises whether these 
steady states can undergo a Turing bifurcation. It is shown in Appendix IbI that 

1. The semi-trivial steady states cannot undergo a Turing bifurcation. 

2. The nontrivial steady state, i?iii, can undergo a Turing bifurcation, but this is so only if 

&23&32 > 1, (7) 

i.e., the mutual inhibition between U2 and U3 must be sufficiently strong — so strong, in particular, 
that they cannot coexist in the absence of the inhibitor (see Figure Et). 

Both conclusions are a consequence of the following fact which will be discussed in Sectional In this model, 
the only destabiHzing mechanism driving the Turing instability is mutual inhibition of U2 and U3, which 
requires the existence of both activators. Thus, the semi-trivial steady states fail to undergo a Turing 
bifurcation because one of the two activators is absent at such steady states. The non-trivial steady state, 
which is characterized by positive concentrations of U2 and U3, allows for a Turing instability, but only if 
their mutual inhibition is sufficiently strong. 
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Although the mutual inhibition must be sufficiently strong to ensure that Em undergoes a Turing 
instability, intuition suggests that it cannot be too strong, lest the two activators become incompatible even 
in the presence of the inhibitor. This is indeed the case. To see this, observe that Q can be satisfied in 
three different ways (see Figure Eb): 

1- &23,^32 > 1, i.e., the mutual inhibition between U2 and U3 is so strong that they display bistable 
dynamics in the absence of the inhibitor. 

2. &23 < 1, ^32 > 1, i.e., in the absence of the inhibitor, U2 inhibits U3 much more than U3 inhibits U2, so 
that U2 ultimately prevails over C/3. 

3. &23 > 1, ^32 < 1, i.e., in the absence of the inhibitor, U3 inhibits U2 much more than U2 inhibits U3, so 
that U3 ultimately prevails over 1/2- 

It turns out that in the first case, when both U2 and U3 inhibit each other strongly, the necessary condition 
Q is satisfied. Yet, the Turing instability cannot be realized because a stable coexistence steady state fails 
to exist even in the presence of the inhibitor. This follows from the fact that Em exists and is stable only 
if the mutual inhibition is sufficiently weak in the presence of the inhibitor, i.e., 

023+121013 , a 032+031012 , 

P23 = ; < -L, P32 = ; < -L 

033 + 031013 022+021012 

which can be recast in the form 

— (&32 - 1) < 021 - 031 < ^ (1 - 623) . (8) 
O12 013 

Clearly, l|Hl cannot be satisfied if 623, 632 > 1- Under these conditions, the mutual inhibition of U2 and U3 is 
so strong that they cannot coexist stably even in the presence of the inhibitor. 

We conclude that Em exists and bifurcates via a Turing instability only if the interaction between U2 
and C/3 in the absence of the inhibitor is such that only one of them prevails ultimately. Furthermore 

1. If [^2 prevails over U3 in the absence of the inhibitor (623 < 1, ^32 > 1), then a Turing instability obtains 
only if 623^32 > 1 and 

< — (632 - 1) < 021 - 031 < — (1 - 623) • (9) 

O12 013 

2. Conversely, if C/3 prevails over U2 in the absence of the inhibitor (623 > 1,632 < 1), then a Turing 
instability obtains only if 623632 > 1 and 

^ (632 - 1) < 021 - 031 < ^ (1 - 623) < 0. (10) 
012 013 

These conditions have a simple physical interpretation. Consider, for instance, the condition (EJ. It says that 
if U2 prevails over C/3 in the absence of the inhibitor, a Turing instability obtains only if (a) 021 — 031 > 0, 
i.e., Ui inhibits U2 more than it inhibits C/3. The stronger inhibition of U2 by C/i is necessary in order to 
compensate for its intrinsic superiority over C/3 . (b) The magnitude of this difference must be neither too 
small nor too large to prevent under- or over-compensation that would preclude the coexistence of U2 and 
C/3 even in the presence of the inhibitor. 

It is shown in Appendix El that the conditions ^ or ljl()|l are not only necessary but almost sufficient 
for Em to undergo a Turing instability. It suffices to impose the additional condition that di be sufficiently 
larger than c?2 and c?3. 

2.4 Simulation of experiments 

To simulate the data shown in Figure QJd, wherein the frontness and backness components segregate spon- 
taneously in response to a uniform chemoattractant profile, we chose parameter values satisfying (jSj, shown 
in the column labeled "C/2 wins" of Table Linear stability analysis shows that given these parameter 
values, the coexistence steady state, Em, undergoes a Turing instability at r ~ 0.6 and wavenumber fc ~ 1 
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Table 1: Parameter values for the simulations. The diffusion coefficient of the inhibitor, Di, is assumed 
to much larger than the diffusion coefficients, 02,0^, of the frontness and backness pathways. The rate 
constants, Ri,Aij, were chosen in order to satisfy the conditions lO or lllOII . These two cases are shown in 
the columns labeled column labeled "J72 wins" and "J/s wins", respectively. 



Parameter 


U2 wins 


U3 wins 


Parameter 


U2 wins 


U3 wins 


Di 


3 


3 


Ai3 


1/3 


4 


D2 


0.001 


0.001 


A21 


3 


1 


Di 


0.005 


0.005 


A22 


1 


1 


Ri 


1 


1 


A23 


1 


2 


R2 


2 


1 


A31 


1 


3 


Rs 


1 


2 


A32 


2 


1 


A12 


4 


1/2 


A33 


1 


1 



(Figure Et)- Computations with the continuation software package CONTENT [2£j confirm the existence 
of this instability (FigureC^). The homogeneous steady state is stable for < r < 0.6, and Turing unstable 
therafter. The concentration profiles of the non-homogeneous steady state created at the Turing bifurca- 
tion point are consistent with the data. Figure 03 shows that the non-homogeneous steady state at r = 1 
is such that U2 and Ui are in phase, and U2 and U3 are out of phase. The latter is consistent with the 
data in Figure We refer to this steady state as the chemoattraction steady state for the following rea- 
son. To a first degree of approximation, the steady state obtained in the presence of a receptor gradient, 
r(x) = 1-1- ecos(7rx),0 < e <g; 1, is identical to the steady state shown in Figure 03. Thus, the profile of 
the frontness pathways is in phase with the distribution of the active receptors, which is characteristic of 
chemoattraction. 

Figures [21 and show the redistribution of the frontness and backness pathways in response to inhibition 
and activation of the frontness pathways. To simulate these experiments, we computed the variation of 
the chemoattraction steady state at r = 1 with respect to i?2. Figure |Ht shows that as R2 decreases, so 
does u^^^ until it becomes zero, i.e., the chemoattraction steady state ceases to exist, and is replaced by 
the homogeneous steady state, -Eioi, in which U2 = throughout the spatial domain. The variation of the 
concentration profile for the frontness pathway is consistent with the data shown in Figure [2| Figures |HJd 
shows that as R2 decreases, so does the gradient of the frontness molecules until it vanishes at R2 ~ 1.5, which 
is in qualitative agreement with the data shown in Figure |2ji. Conversely, at sufficiently large values of R2, 
the non-homogeneous steady state is replaced by the homogeneous steady state, Euq, characterized by high 
levels of the frontness molecules uniformly distributed throughout the spatial domain (see Figure Eb) . The 
variation of the concentration profile for the backness pathway is also consistent with the data. Figures |HJ; 
shows that as R2 decreases, the backness pathway progressively advances into the pre-existing front until it 
occupies not only the back but also the front of the cell. This is consistent with the data shown in FigureEt- 

Simulations in which r is held at 1.0 and R3 is progressively changed give analogous results (simulations 
not shown). As R3 increases, the frontness pathways occupy a progressively smaller portion of the cell until 
there is no frontness throughout the cell. The simulations did not yield multiple peaks corresponding to the 
multiple pseudopods shown in Figure^. Numerical calculations of the critical wavenumber, fco, show that 
it increases as R3 decreases, which is the correct trend. However, the steady state ceases to exist before ko 
can exceed 2. 

To investigate the existence of steady states corresponding to chemorepulsion, we chose parameter values 
satisfying IjlOfl so that U3 prevails over U2 in the absence of the inhibitor (column labeled "J/s wins" of 
Table nj. Here, the coexistence steady state, Em, undergoes a Turing bifurcation at r « 0.65 and k — 1 
(Figures Eb)- The non-homogeneous steady state that emerges from this bifurcation is such that U2 and 
U3 are once again out of phase, ^ but it is U3, instead of U2, that is in phase with Ui (Figure Et) . This 
non-homogeneous steady state represents the chemorepulsion steady state, since it approximates the steady 

^It is shown in Appendix ^ that regardless of the parameter values, U2 and U3 will always be out-of-phase in the non- 
homogeneous steady state bifurcating from Em. 
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Figure 7: Symmetry-breaking of the non-trivial steady state when the parameters are chosen such that 
U2 wins in the absence of the inhibitor, (a) The bifurcation diagram depicting the variation of u™^^ = 
maxo<a;<i U2(x) as a function of the active receptor concentration, r. The thin hne represents the homoge- 
neous non-trivial steady state, -Em; the dashed portion of this line indicates that Em is Turing unstable. 
The thick line represents the non-homogeneous steady state emerging from the Turing bifurcation at r « 0.6. 
(b) The profiles of ui, U2, and U3 at the non-homogeneous steady state corresponding to r = 1. The steady 
state was attained by perturbing the homogeneous steady state at r = 1, i.e., by integrating equations I^HSj 
with initial conditions, ui(Q,x) = ui, Ui{0,x) = Ui[\ + 0.01 cos(7rx)], i — 2,3, where Ui = 1,2,3 denotes the 
non-trivial homogeneous steady state at r = 1. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 oL . . , . I . , , , , 

R 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) (b) (c) 



Figure 8: Variation of the chemoattraction steady state with respect to i?2 at r = 1. (a) The bifurcation 
diagram showing the variation of u™^^ as a function of i?2 • The thick line represents the non-homogeneous 
steady state. The thin line represents the homogeneous semi-trivial steady state, i^ioi- As R2 decreases, so 
does u™^^ until it intersects the branch, uf^^ = 0, which corresponds the homogeneous semi-trivial steady 
state, Eioi. (b) As R2 decreases, the concentration profile of U2 decreases in magnitude and occupies a 
progressively smaller proportion of the cell until U2 is identically zero, (c) As R2 decreases, the concentration 
profile of increases in magnitude and advances into the pre-existing front until U3 acquires the uniform 
profile of i?ioi. 
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Figure 9: Symmetry-breaking of the non-trivial steady state when the parameters are chosen such that 
C/3 wins in the absence of the inhibitor, (a) The bifurcation diagram depicting the variation of u™*^ = 
maxo<i,<i 1/3(2;) as a function of the active receptor concentration, r. The thin Une represents the ho- 
mogeneous non-trivial steady state, Em; the dashed portion of this line indicates that Em is Turing 
unstable. The thick line represents the non-homogeneous steady state emerging from the Turing bifurcation 
at r « 0.65. (b) The profiles of ui, U2, and M3 at the non-homogeneous steady state attained by perturbing 
the homogeneous steady state at r = 1. 



state that would be obtained in the presence of a receptor gradient, r{x) = 1 + ecos(7ra;),0 < e <C 1, and 
the profile of the frontness pathways are, in this case, out of phase with the active receptor distribution. 
Given the symmetry of the equations, it is not surprising that the model admits a chemorepulsion steady 
state. For, if choose the parameters corresponding to the column wins" in Tabled but interchange U2 
and M3, the chemoattraction steady state in Figure 03 turns into a chemorepulsion steady state. However, 
the analysis of the model suggests a mechanism for the transition. Specifically, the conditions l|9 HT(Hl imply 
that transitions from chemorepulsion to chemattraction, such as those observed in neurons [22|, can be 
triggered by pharmacological agents that selectively inhibit the backness pathway so that the frontness 
pathway becomes intrinsically superior in the absence of the inhibitor. 



3 Discussion 

We have shown that a model based on mutual inhibition of the frontness and backness pathways can yield 
spontaneous polarization involving spatial segregation of the two pathways. The frontness and backness 
pathways are mutually incompatible in the absence of the inhibitor (eq. , but coexist in the presence of 
the inhibitor because it suppresses the growth, and hence, the mutual incompatibility, of the frontness and 
backness pathways. Since the inhibitor is diffusible, it can suppress the nominal incompatibility manifested 
at low receptor levels. However, it fails to achieve this goal at high receptor levels, resulting in spatial 
segregation of the frontness and backness pathways (see Figures and . 

Since the model takes due account of the backness pathways, we could explore the effect of frontness 
inhibition or activation on the backness pathways, and the effect of backness modification on the frontness 
and backness pathways. The simulations in Section [21 yielded results that are in qualitative agreement with 
the data shown in FiguresQ}^ These results are beyond the scope of activator-inhibitor or activator-substrate 
models of gradient sensing since they do not account for the backness pathways. 

The particular manner in which the two pathways segregate depends on their behavior in the absence of 
the tempering infiuence of the inhibitor. If the frontness pathways dominate over the backness pathways, 
their spatial segregation is consistent with chemoattraction. Conversely, if the backness pathways dominate. 
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the emergent spatial pattern corresponds to chemorepulsion. The model therefore suggests that transitions 
between chemoattraction and chemorepulsion can be provoked by agents that shift the balance of power 
between the frontness and backness pathways. 

The model proposed here is the simplest possible model of spontaneous polarization driven by mutual 
inhibition. By this, we mean the following 

1. A model containing a smaller number of variables will not yield spontaneous polarization. 

2. The model is simpler than any other model in which the symmetry-breaking is driven entirely by 
mutual inhibition. 

In what follows, we justify these conclusions, since they are required to explain certain points raised above. 

To see that a model containing 2 variables will not yield spontaneous polarization driven by mutual 
inhibition, let u = 1*2]* be a homogeneous steady state of a two-component reaction-diffusion system 

du, r, du, 

'dt =^'9;;2 '"2), ^= 1,2 

and let J{u) denote the Jacobian at u. Then, u can undergo a Turing instability only if jJI] 

Jii + J22 < 0, Jii J22 - J12 J21 > 0, D1J22 + D2J11 > 0. (11) 

The first and third conditions imply that Jn J22 < 0, i.e., one of the components, say Ui, must inhibit its 
own synthesis (Jn < 0), and the other, say U2, must activate its own synthesis (J22 > 0). Now, if the 
two components inhibit each other, then J12, J21 < 0, so that J12J21 > and Jn J22 — J12J21 < 0, which 
violates the second condition for Turing stability. Thus, a two-component model in which each component 
inhibits the other component cannot yield spontaneous polarization. Indeed, such a model can yield a Turing 
instability only if J12J21 < 0, i.e., the interaction between Ui and U2 is asymmetric. If Ui inhibits U2, then 
U2 must activate Ui; this corresponds to the activator-inhibitor model; conversely, if Ui activates U2, then 
U2 must inhibit Ui, which corresponds to the activator-substrate model. No other cross-interactions are 
admissible. It follows that if spontaneous polarization is to be driven by mutual inhibition, there must be 
at least three components. Thus, we were led to modify the Bourne model by postulating the existence of 
an additional component, namely, a diffusible inhibitor. 

To show that the model is simpler than any other model in which symmetry-breaking is driven entirely by 
mutual inhibition, we appeal to the fact that in a general reaction-diffusion system, a homogeneous steady 
state cannot become Turing instable unless it contains an unstable subsytem [2HI-^ The three-component 
model considered here has three 1-component subsystems, {C/i, C/2, U^}, and three 2-component subsystems, 
{UiU2,U2U^,UiUj} . It turns out that in the neighborhood of the non-trivial steady state, Em, all three 
1-component subsystems are stable, i.e. Jn, J22, J33 < fsee lT^jl . In other words, the self-interactions of the 
individual components are stabilizing. Notably, U2 and are stable subsystems even though their synthesis 
was assumed to be autocatalytic. This is because the kinetics of their synthesis are such that the feedback is 
positive only when U2 and M3 are small. However, in the neighborhood of -Em, U2 and M3 are so large that the 
self-limiting effect dominates, and the feedback is, in fact, negative. Now, since the 1-component subsystems 
are stable, Em can be Turing unstable only if at least one of the 2-component subsystems is unstable, i.e., 
the cross-interaction of at least one 2-component system is destabilizing. Such destabiHzing cross-interactions 
can occur only if the interaction between the 2 components is mutually inhibitory or synergistic, i.e., 

J12J21 > or J13J31 > or J23 J32 > 0. 

The first two conditions cannot be satisfied since the interaction between Ui and each of the two activators, 
U2 and is neither mutually inhibitory nor mutually synergistic — Ui inhibits U2 and U^, but U2 and 

subsystem refers to any proper subset of the model variables. For example, the two-component model discussed above 
has two subsystems, {Ui, U2\- A subsystem is said to be unstable if it becomes unstable as soon as the remaining components of 
the system are somehow rendered constant. In other words, without the stabilizing effects of the other variables, the subsystem 
diverges from its values at the homogeneous steady state. In the 2-component model, the activator is an unstable subsystem: 
If the inhibitor or substrate level is somehow held constant, the activator diverges from the homogeneous steady state due to 
positive feedback (J22 > 0). As shown above, the existence of such an unstable subsystem is necessary for the full system to 
undergo a Turing instability. 
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activate Ui. Thus, {Ui, U2} and {Ui, U3} are stable subsystems, and {U2, U3} is the only unstable subsystem 
in the model. Consequently, the spatial segregation of the frontness and backness pathways is driven entirely 
by mutual inhibition between U2 and C/3 — positive feedback plays no role since Jn, J22, J33 < 0. 

It should be noted that reciprocal distribution of the frontness and backness pathways (specifically, PI3K 
and PTEN) is a central feature of the local-excitation-global-inhibition model proposed by Iglesias and 
coworkers However, the mechanism for spatial segregation of PI3K and PTEN is such that polarization 
occurs only in the presence of a chemoattractant gradient. Indeed, according to this model 

1. Exposure to chemoattractant generates binding sites for PI3K and PTEN at the membrane. 

2. The rate of generation of binding sites for PI3K and PTEN is proportional to and r, respectively, 
where r denotes the active receptor concentration. 

3. The concentration of PI3P is determined by the relative concentrations of membrane-bound PI3K and 
PTEN — the higher the PI3K:PTEN ratio, the larger the concentration of PI3P. 

It follows from (2) that when a cell is exposed to a gradient, the leading edge develops a high PI3K:PTEN 
ratio, while the trailing edge has a high PTEN:PI3K ratio. Then, (3) implies that the PI3P concentration 
at the leading edge is higher than that at the trailing edge. Importantly, the polarized distribution of 
PI3P arises only in the presence of a gradient. For, if r is constant, the PI3K:PTEN ratio, and hence, the 
concentration of PI3P, is constant all over the cell membrane. Thus, the local-excitation-global-inhibition 
model cannot capture the spontaneous polarization observed in the absence of gradients. 

Although our model captures spontaneous polarization driven by mutual inhibition of frontness and 
backness pathways, it is missing two important features of the experimental data. First, it is known that 
both negative and positive feedback operate in chemoattractant-mediated polarization (see Figure and 
refs. [2310]). Although we assumed the existence of positive feedback in the model, it played no role in 
the spontaneous polarization since J22, J33 < in the neighborhood of Em. Second, the model does not 
account for adapation, the process that enables the cell to ultimately return to the homogeneous steady state 
even in the presence of chemoattractant [221 • A modified model taking due account of positive feedback and 
long-term adaptation is currently under investigation. 

Two of the model variables {U2, U3) correspond to the frontness and backness pathways of the Bourne 
scheme. However, the model hypothesizes the existence of an additional variable, namely, a diffusible in- 
hibitor required to suppress the mutual incompatibility of the frontness and backness pathways. It is, by 
no means, necessary that this diffusible component be an inhibitor. It could just as well be a diffusible 
substrate that is required for the synthesis of the two activators, but is consumed in process of activator 
synthesis. This would reverse the signs of J12, J13, J21, J31, but the key inequalities in the above arguments, 
J13J31, J12J21 < 0, would be preserved. A diffusible inhibitor was assumed only because there is some 
evidence of their existence. 

1. Luo and coworkers have shown that cytosolic inositol phosphates, whose levels surge immediately after 
chemoattractant stimulation, inhibit the frontness pathways in both Dictyostelium (201 neutrophil- 
like HL-60 cells [S]. 

2. Nimnual et al have shown that Rac-mediated production of reactive oxygen species downregulates the 
backness component, Rho [32 . 

It remains to be seen if these diffusible inhibitors play a critical role in the maintenance of polarity, as 
required by the diffusible inhibitor of the proposed model. 

4 Conclusions 

We have formulated a model that provides a mathematical realization of the Bourne scheme for sponta- 
neous polarization by mutual inhibition of frontness and backness pathways [15 • The model predicts several 
experimentally observed features that are outside the scope of prevailing models of gradient sensing. 
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1. Mutual inhibition of the frontness and backness pathways plays a critical role in generating the spon- 
taneous polarization. 

2. If the frontness pathway is suppressed, the backness pathways occupy a progressively larger proportion 
of the cell. Conversely, if the backness pathway is suppressed, the frontness pathways invade and 
occupy the previous back of the ceh. 

3. Depending on the parameter values, the model yields steady states corresponding to both chemoat- 
traction and chemorepulsion. Analysis of the model suggests that chemorepulsion-to-chemoattraction 
transitions observed in neurons can be triggered by agents that inhibit the backness pathways more 
than the frontness pathways. 

The model provides a useful starting point for formulating models that account for the positive and negative 
feedback effects, both of which have been shown to play a role in gradient sensing. 
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A Existence/stability of homogeneous steady states 

The homogeneous steady states of (0HSl satisfy the equations 

—ui + ai2U2 + ai3U3 = 
(r - a2iui ~ a22U2 — 023^3) U2 — 
(r - a3iui - a32U2 ~ 033^3) U3 = 0. 

It follows that there are 4 possible steady states: ui = U2 = U3 = 0, denoted i?ooo; ui,U2 > 0,U3 — 0, 
denoted -Buo; "1,^3 > 0,U2 = 0, denoted i?ioi; and 1*1,^2,^3 > 0, denoted Em. In what follows, we show 
the existence and stability criteria for all 4 steady states. To this end, it is convenient to note that the 
Jacobian at any point, ui,U2,it3 is 



J = 



-1 

-p2U2a21 



ai2 

-p2U2a22 + P2/2 
-P3"3a32 



ai3 

-/92W2a3i 
-p3U3a33 + P3/3 



where 



/2 
/3 



r — a2iui 
r - a3iui 



022^2 - a23U3 
032^2 - a33W3 



1. i?ooo = [0,0,0]*: It is evident that this steady state always exists. The Jacobian yields no information 
regarding the stability since it is singular at Eqoo- However, we can infer that it is unstable for all 
r > by observing that the U2-axis is an invariant manifold. The motion along this invariant manifold 
is given by 

du2 I . 

—— = p2[r - a22U2)U2, 

at 

which drives the system away from Eqoq for all r > 0. 

2. i?iio = [ui,U2, 0]', ui,U2 > 0: It is easy to check that this steady state exists for all r > 0, and is given 

by 

r 

Us = 0, U2 = , Ul = ai2U2, 

a22 



where 022 = 022 + 021012, and the Jacobian at -Eiio is 



-1 ai2 

-P2U2a21 -p2U2a22 





It follows that one of the eigenvalues is 



ai3 
-p2U2a3i 



a32 



(12) 



Ai = psir - asiui - 032^2) = Psr 1 - 



"33 



where Q!32 = 032 + 031012. The other two eigenvalues satisfy 

A2 + A3 = -1 - p2a22ru2 < 0, 

A2A3 = p2a22ru2 + P2U2a2iai2 > 0, 

so that the real parts of A2, A3 are negative. Hence, i?iio is stable if and only if /332 

3. i?ioi ~ [uijO, U3]*, ui,U3 > 0: The calculations are analogous to those for i?ioi. 

4. Em ~ 1*2,^3]*, Mi,U2,M3 > 0: This steady state satisfies the equations 



a32/a22 > 1- 



Ul 

r 
r 



ai2U2 + ai3it3 

a22U2 + Ce23U3 
a32U2 + Q;33M3. 
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It follows that Em exists for all r > if and only if [1, 1]* lies between [022, 032]* and [023, 0:33]*, i.e., 
/?23 = a23/a33,^32 = 032/022 Satisfy 

/323,/^32 < 1 or /323,/?32 > 1- 

In both cases, the steady state is given by 

1 1-/323 

U2 = r: r, r, 

Q!22 -L — P23P32 

1 I-P32 
O33 -L — P23P32 

ui = ai2U2 + aisua. 



Evidently, ui,M2,U3 increase linearly with r. 

It turns out that Em is unstable if /323,/932 > 1, and stable for sufficiently small r > if /?23,/332 > 1- 
To see this, observe that the Jacobian at Em is 



-1 ai2 ai3 

-^2^2021 -^2^2022 -p2U2a31 
-^3^3031 -^3^3032 -p3U3a33 



(13) 



so that 



tr J 

detJ 



-1 - p2a22U2 - P3a33W3, 
-P2P3W2W3a22a23(l " /323/332), 
S J = P2P3U2U3{a22a33 - 023032) + ^2^2022 
+ P3U30:33- 

If /323,/332 > 1, then det J > 0, and Em is unstable. Hence, Em is stable only if /323,/332 < 1- We 
are particularly interested in the case when 022033 — 023O32 > 0, since this is necessary for Turing 
instability (see below). Under these conditions, Em is stable for all sufficiently small r > 0, since 
tr J < 0, det J < for all r > 0, and (tr J)(S J) - det J < for all sufficiently small r > fFigurelTnll. 



B Turing instability of homogeneous steady states 
B.l Conditions for Turing instability 

A homogeneous steady state, u = [ui{r),U2{r),U3{r)Y , is Turing unstable if it is stable in the absence of 
diffusion and unstable in the presence of diffusion. Now the stability of the homogeneous steady state, u, in 
the presence of diffusion is determined by the linearized equation 

Ov Dv 

- = n^^JM. ,14) 

where v denotes the deviation from u, J{u) denotes the Jacobian at u, and D = diag(c?i, da, dg). Let fik and 
4>k{x) denote the eigenvalues and eigenfunctions of the differential operator, with Neumann boundary 
conditions, i.e., 

P-k — —{kir)'^, 4'k{x) — cos(fc7ra::), k — 0,1,2, . . . 

If the steady state is perturbed along any eigenfunction, ipkix), the subsequent evolution of the perturbation 
is of the form v{t, x) — Sk{t)(f>k{x) where Sfc(t) G R^. Substituting in ifnt yields 

^ = CkSk, Ck{u) = J{u) - {kirfD. 
at 

It follows that the homogeneous steady state is stable if and only if the eigenvalues of Ck{u) have negative 
real parts for all A; > 0. It is Turing unstable if and only if 
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Figure 10: Variation of tr J ( — ), det J ( — ), and (tr J)(SJ) — det J ( ) as a function of r. (a) Parameter 

values chosen such that Ui prevails over U-i in the absence of U\. (b) Parameter values chosen such that 
C/3 prevails over U2 in the absence of J7i. 

1. The eigenvalues of Co(u) — J{u) have negative real parts (which ensures that the homogeneous steady 
state is stable in the absence of diffusion) . 

2. The eigenvalues of Ck{u) have positive real part for some A: > (which ensures that the homogeneous 
steady state is unstable in the presence of diffusion). 

The Turing bifurcation occurs at the critical value, r = ro > 0, such that exactly one of the eigenvalues of 
Ck{u) becomes zero, i.e., 

d 

for some fc > (Figure EJ- 



det Cfc(u) = 0, -j^ det Ck{u) = 



B.2 Turing instability of the homogeneous steady states 

It follows from the definition of Ck that 

det Ck{u) = Ai23 - (diA23 + d2Ai3 + d3Ai3)k^ 

+ (c?id2A3 + d2d3Ai + di(i3A2)fc4 

- (dirf2rf3)fc' (15) 

where A^, Ay,z ^ j, A123 are the determinants of the 1-, 2-, and 3-dimensional subsystems of J, i.e., 

Aj = Jji(u), 



Aij = det 



Jii Jij 
Jji Jjj 



and A123 — det J(m). These results are sufficient for investigating the possibility that the homogeneous 
steady states undergo a Turing bifurcation. 

We consider each of the homogeneous steady states. 

1. -Eooo = [0,0,0]*: Since i?ooo is unstable even in the absence of diffusion, it cannot entertain Turing 
instability. 
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Figure 11: Variation of det Ck as a function of at various values of r. (a) Parameter values chosen such 
that U2 prevails over U3 in the absence of Ui. (b) Parameter values chosen such that C/3 prevails over U2 
in the absence of J7i. 



2. Eiio = [mi,M2,0]*: This steady state is stable in the absence of diffusion provided (3^2 > 1- Thus, it 
can undergo a Turing instability if and only if there is an r > such that ifT^ll has a positive root. It 
turns out that such a root does not exist. To see this, observe that since Euq is stable in the absence of 
diffusion, A123 = det J < 0. Inspection of the Jacobian II12|I shows that the A.^'s are always negative, 
and the A^ 's are always positive. Hence, all the coefficients of the polynomial Ijl5|l are negative for all 
r > 0. Such a polynomial cannot have positive roots (fc > 0). We conclude that Euq cannot undergo 
a Turing instability. 

3. -Eioi = [ui, 0,1*3]*: This steady state is stable in the absence of diffusion provided P23 < 1- An 
argument analogous to that for Euq shows that it cannot undergo a Turing instability. 

4. Em — [ui, 1*2,1*3]*: This steady state is stable in the absence of diffusion provided /?23,/?32 < 1 and 
r is sufficiently small. Under these conditions, A123 < 0. Inspection of the Jacobian Ijl3|l shows that 
no matter what the value of r, all the A^'s are negative, and A12, A13 are positive. Hence, Em can 
undergo a Turing instability only if A23 < 0. This condition is not only necessary, but also sufficient. 
For, if d2 and c?3 are sufficiently small compared to di, the coefficient of is approximately rfi A23 < 0, 
and the polynomial lfT5|l has a positive root. 

The non-homogeneous steady state that bifurcates from Em is such that U2 and C/3 are out-of-phase. 
To see this, it suffices to observe that to a first degree of approximation, the non-homogeneous steady 
state near the Turing bifurcation point has the form vq cos{kQx), where kg be the critical wavenumber, 
and Wo is the eigenvector corresponding to the zero eigenvalue of Cko, i-e., CkoVo = 0- The null space 
of the rank-2 matrix, Cko, is spanned by the vector, [ci, C2, 1], where 

_ Cko,22 — C ka ^12C ko, 21 /C kg m 
^2 ^ I 

L-fco,23 — L'feo,21<-^feo,13/<-'feo41 

Substitution of the coefficients of Cko shows that C2 is always negative, which implies that C/2 and C/3 
are out-of-phase. 
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